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i  A.-  Statesoent  of  Work 

The  focus  of  this  research  project  has  been  the  study  of  systems  of 
conservation  laws  that  lose  strict  hyperholicity.  The  approach  has  been  to 
consider  two  areas  of  application#  namely  isentropic  gas  dynamics  for  a  van 
der  Waals  gas#  and  multiphase  flow  in  a  porous  medium.  Additionally#  a 
numerical  study  of  the  elastic  string  equations  using  Gliinn's  random  choice 
method  has  been  completed. 

I 

B.  Summary  of  Results 

1.  I  have  considered  the  role  of  undercompressive  shocks#  or  phase 
jumps#  insolving  the  initial  value  problem  for  isentropic  gas  dynamics 
with  a  nonmonotone  equations  of  state.  The  main  results  are: 

(a)  If  the  initial  data  is  close  to  the  Maxwell  line  of  equilibrium 
thermodynamics#  then  the  Riematvi  initial  value  problem  has  a 
centered  solution  consisting  of  a  fast  wave  and  a  slow  wave 
separated  by  a  slow  phase  jump  (see  [1]). 

(b)  For  certain  initial  data#  the  Riemann  problem  has  two 
solutions#  distinguished  by  the  presence  or  absence  of  a  pair 
of  phase  jumps.  This  is  summarized  in  [2]. 

2.  The  equations  of  motion  of  three  phase  flow  in  a  porous  medium 
constitute  a  2x2  system  of  conservation  laws.  In  a  recent  paper  [3]#  I 
used  a  topological  argument  to  show  that  these  equations  are  generally  of 
mixed  elliptic  and  hyperbolic  type.  The  elliptic  part  of  state  space  is 
a  bounded  region. 

The  main  part  of  this  project  has  been  the  study  of  equations  with 
isolated  umbilic  points.  These  correspond  to  elliptic  regions  shrunk  to 
a  point.  In  [4]#  we  classified  equations  with  umbilic  points  by 
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obtaining  a  normal  form  for  equations  with  quadratic  nonlinearities.  The 
equations  of  porous  media  flow  fall  into  two  of  the  four  classes/  whereas 
all  four  classes  have  been  found  in  applications  in  elasticity/  in  which 
there  is  additional  symnetry.  In  [5]/  we  solved  the  Riemann  problem  for 
one  case  having  symnetry/  introducing  undercompressive  shocks/  cmd  idee^ 
from  bifurcation  theory.  In  [6]/  the  bifurcation  ideas  are  developed 
further  to  solve  the  Riemann  problem  in  three  of  the  four  Ccises.  This 
paper  also  introduces  and  classifies  boundaries  of  regions  in  state  space 
across  v^ich  the  solution  of  the  Rienann  problem  undergoes  a  qualitative 
cheuige  such  as  a  chauige  in  the  type  of  waves  in  the  solutions. 

One  objective  of  this  research  has  been  to  understand  the  solution 
of  the  Riemann  problem  sufficiently  well  that  an  algorithm  can  be 
developed  to  automate  the  solution.  lliis  is  now  being  implemented  by 
conputing  the  wave  curves  (themselves  boundaries)/  which  undergo 
transitions  across  boundaries.  The  classification  of  boundaries  plays  an 
important  role  in  the  computations. 

3.  The  implementation  of  Glimm's  scheme  for  the  4x4  syst«n  of  first 
order  equations  describing  the  motion  of  an  elastic  string  has  been 
carried  out  in  [7].  The  algorithm  uses  the  solution  of  the  Riemann 
problem  that  I  had  worked  out  earlier  in  a  form  that  is  computationally 
viable.  The  boundary  conditions  are  handled  by  solving  the  Goursat 
(quarter  plane)  problem/  which  has  a  subtlety  due  to  the  possibility  of 
the  string  doubling  over.  The  most  significant  part  of  this  work  is  the 
observed  link  between  new  explicit  solutions  of  the  equations/  involving 
purely  transverse  motion/  and  a  numerical  solution  that  appears  to  be 
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periodic*  and  almost  entirely  transverse.  This  is  significant  because  it 
is  generally  thought  that  shocks  form  from  spatially  periodic  initial 
data*  but  no  theorems  are  known  for  4x4  systems. 
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